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Using nuclear magnetic resonance techniques, we experi- 
mentally investigated the effects of applying a two bit phase 
error detection code to preserve quantum information in nu- 
clear spin systems. Input states were stored with and with- 
out coding, and the resulting output states were compared 
with the originals and with each other. The theoretically ex- 
pected result, net reduction of distortion and conditional er- 
ror probabilities to second order, was indeed observed, despite 
imperfect coding operations which increased the error proba- 
bilities by approximately 5%. Systematic study of the devia- 
tions from the ideal behavior provided quantitative measures 
of different sources of error, and good agreement was found 
with a numerical model. Theoretical questions in quantum er- 
ror correction in bulk nuclear spin systems including fidelity 
measures, signal strength and syndrome measurements are 
discussed. 



I. INTRODUCTION 

Recent progress in experimental implementation of 
quantum algorithms has demonstrated in principle that 
quantum computers could solve specific problems in 
fewer steps than any classical machine These first 

generation quantum computers were 2-spin molecules in 
solution. They were initialized, manipulated and mea- 
sured at room-temperature using bulk nuclear magnetic 
resonance (NMR) spectroscopy techniques Clas- 
sical redundancy in the ensemble and the discrete na- 
ture of the answers ensured that the correct answers 
were obtained despite gate imperfections and moderate 
rates of decoherence. However, the accumulation of er- 
rors would be detrimental in larger quantum computers 
and for longer computations in the future, and methods 
to protect quantum information will be needed. 

Classically, errors (bit flips) are detected and fixed by 
error correcting codes. Information is encoded redun- 
dantly, and the output contains information on both the 



encoded input and the errors that have occurred, so that 
the errors can be reversed. Generalization to quantum 
information is non-trivial since it is impossible to clone 
an arbitrary quantum bit (qubit) and to measure quan- 
tum states without disturbing the system. Furthermore, 
there is a continuum of possible errors, and finally, entan- 
glement can cause errors to propagate rapidly throughout 
the system. 

Despite the apparent difficulties, quantum error cor- 
rection was shown to be possible theoretically, and can 
be useful for reliable computation even when coding op- 
erations are imperfect. Shor and Steane pO|,pT[ realized 
a major breakthrough by constructing the first quantum 
error correcting codes. Prudent use of quantum entan- 
glement enables the information on the errors to be ob- 
tained by non-demolition measurements without disturb- 
ing the encoded inputs; it also enables digitization of er- 
rors. These schemes correct for storage errors but not 
for the extra errors introduced by the coding operations. 
The extension to handle gate errors and to achieve reli- 
able computation with a certain accuracy threshold were 
subsequently developed by many others |]l2|-|l7||. 

In this paper, we report experimental progress toward 
this elusive goal of continued quantum computation. Wc 
implement a simple phase error detection scheme |lq] 
that encodes one qubit in two and detects a single phase 
error in either one of the two qubits. The output state 
is rejected if an error is detected so that the probability 
to accept an erroneous state is reduced to the smaller 
probability of having multiple errors. Our aim is to 
study the effectiveness of quantum error correction in a 
real experimental system, focusing on effects arising from 
imperfections of the logic gates. Therefore, the experi- 
ment is designed to eliminate potential artificial origins 
of bias in the following ways. First, we compare out- 
put states stored with and without coding (the latter 
is unprotected but less affected by gate imperfections). 
Second, by ensuring that all qubits used in the code de- 
cohere at nearly the same rate, we eliminate apparent 
improvements brought by having an ancilla with lifetime 
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much longer than the original unencoded qubit. Third, 
our experiment utilizes only naturally occuring error pro- 
cesses. Finally, the main error processes and their rela- 
tive importance to the experiment are thoroughly stud- 
ied and simulated to substantiate any conclusions. In 
these aspects, our study differs significantly from previ- 
ous work fll9| demonstrating quantum error correction 
working only in principle. 

Using nuclear spins as qubits, we performed two sets 
of experiments: (1) The "coding experiments" in which 
input states were encoded, stored and decoded, and (2) 
the "control experiments" in which encoding and decod- 
ing were omitted. Comparing the output states obtained 
from the coding and the control experiments, both error 
correction by coding and extra errors caused by the im- 
perfect coding operations were taken into account when 
evaluating the actual advantage of coding. In our ex- 
periments, coding reduced the net error probabilities to 
second order as predicted, but at the cost of small ad- 
ditional errors which decreased with the original error 
probabilities. We identified the major imperfection in the 
logic gates to be the inhomogeneity in the radio frequency 
(RF) field used for single spin rotations. Simulation re- 
sults including both phase damping and RF field inhomo- 
geneity confirmed that the additional errors were mostly 
caused by RF field inhomogeneity. The causes and ef- 
fects of other deviations from theory were also studied. 
Tomography experiments giving the full density matrices 
at major stages of the experiments further confirmed the 
agreement between theoretical expectations and the ac- 
tual results. 

The rest of the paper is structured into five sections: 
Section [n] consists of a comprehensive review of the back- 
ground material for the subsequent sections of the paper. 
This background material includes the phase damping 
model, the two bit coding scheme, and the theory of 
bulk NMR quantum co mputa tion. These are reviewed 



in Sections [I A, II B and II C Readers who are familiar 



with these subje cts can skip the appropriate parts of the 
review. Section III describes the methods to implement 
the two bit coding scheme in NMR, and the fidelity mea- 
sures to evaluate the scheme. Section presents the 
experimental details. Section ^ consists of the experi- 
mental results together with a thorough analysis. The 
effects of coding, gate imperfections and the causes and 
effects of other discrepancies are studied in detail. In 
Section VI, we conclude with a summary of our results. 
We also discuss syndrome measurement in bulk NMR, 
the equivalence between logical labeling and coding, the 
applicability of the two bit detection code as a correction 
code exploiting classical redundancy in the bulk sample 
and the signal strength issue in error correction in bulk 
NMR. Sections III and ^ contain the main results of the 
paper; the remainder is included for the sake of complete- 
ness and to develop notation and terminology we believe 
will be helpful to the general reader. 



II. THEORY 

A. Phase damping 

Phase damping is a decoherence process that results 
in the loss of coherence between different basis states. It 
can be caused by random phase shifts of the system due 
to its interaction with the environment. For example, let 
\ip) = a\0) + b\l) be an arbitrary pure initial state. A 
phase shift, P, can be represented as a rotation about 
the z axis by some angle 9, 



P = exp 
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where a z is a Pauli matrix. The resulting state is given 
by P\ijj) = ae- l9 ' 2 \Q) + be ie ' 2 \l). Let p be the density 
matrix of the initial qubit, 
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After the phase shift given by Eq.(Q), the density matrix 
becomes 
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Here, we model phase randomization as a stochastic 
Markov process with 9 drawn from a normal distribution. 
The density matrix resulting from averaging over 9 is 
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where s 2 is the variance of the distribution of 9. By 
Markovity, the total phase shift during a time period t is 
a random walk process with variance proportional to t. 
Therefore, we replace s 2 /2 by Xt in Eq.(Q) when the time 
elapsed is t. Since the diagonal and the off-diagonal el- 
ements represent the populations of the basis states and 
the quantum coherence between them, the exponential 
decay of the off-diagonal elements caused by phase damp- 
ing signifies the loss of coherence without net change of 
quanta. 

Phase damping affects a mixed initial state similarly: 
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since the density matrix of a mixed state is a weighted 
average of the constituent pure states. 

One can represent an arbitrary density matrix for one 
qubit as a Bloch vector (x,y,z), defined to be the real 
coefficients in the Pauli matrix decomposition 
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where the Pauli matrices are given by 
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The space of all possible Bloch vectors is the unit sphere 
known as the Bloch sphere. Phase damping describes the 
axisymmetric exponential decay of the x and y compo- 
nents of any Bloch vector, as depicted in Fig. [|. 
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FIG. 1. Trajectories of different points on the Bloch sphere 
under the effect of phase damping. Points move along perpen- 
diculars to the z-axis at rates proportional to the distances 
to the i-axis. As a result, the Bloch sphere turns into an 
ellipsoid. 



of £(p) are not generally obtainable because they are not 
necessarily orthogonal to each other. 

We emphasize that Eqs.(|^) and (^) describe the same 
physical process. Eq.(^) provides a discrete interpreta- 
tion of phase damping, with the continuously changing 
parameter e~ xt embedded in the probabilities of the pos- 
sible events. 

For a system of multiple qubits, we assume indepen- 
dent decoherence on each qubit. For example, for two 
qubits A and B with error probabilities p a and pb, the 
joint process is given by 

£{ P ) = (1 - p„)(l - p b ) (J <g> /) p (I ® J) 
+ {l-p a )pb (I ® a z ) p (I ® a z ) 

+ p a (1 - p b ) (<T Z 2)1) p (<J Z <g> /) 

+ PaPb (cr z (g) a z ) p (<J Z (g) <J Z ) , (10) 

where p now denotes the 4x4 density matrix for the 
two qubits. The events a z ® I and I ® a z are first order 
errors, while a z <g> a z is second order. First and second 
order events occur with probabilities linear and quadratic 
in the small error probabilities. 

Having described the noise process, we now proceed to 
describe a coding scheme that will correct for it. 

B. The two bit phase damping detection code 



In contrast to the above picture of phase damping as a 
continuous process, we now describe an alternative model 
of phase damping as a discrete process. This will facil- 
itate understanding of quantum error correction. The 
essence is that any quantum process p — > £(p) can be 
written in the operator sum representation as 



£{p) = Y J A k pA\ 



(8) 



where the sum is over a finite number of discrete 
events, A k , which are analogous to quantum jumps, and 
^ k A\A k is positive. For instance, Eq.(||) describing 
phase damping can be rewritten as 



£{p) = (1 -p) IpP + p a z pa 
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wherep= (1 — e~ At )/2. In Eq.(^J), the output £(p) can be 
considered as a (l-p):p mixture of p and a z pa\\ in other 
words, £(p) is a mixture of the states after the event "no 
jump" (J) or "a phase error" (er 2 ) has occurred. The 
weights 1—p and p are the probabilities of the two possi- 
ble events. In general, each term A k pA k in Eq.(|J) repre- 
sents the resulting state after the event A k has occurred, 
with probability tr [AkpA\). This important interpreta- 
tion is used throughout the paper. Note that the decom- 
position of £ (p) is a mathematical interpretation rather 
than a physical process. The component states A k pA k 



Quantum error correction is similar to its classical 
analogue in many aspects. First, the input is encoded 
in a larger system which goes through the decoherence 
process, such as transmission through a noisy channel 
or storage in a noisy environment. The encoded states 
(codewords) are chosen such that information on the un- 
desired changes (error syndromes) can be obtained in the 
extra degrees of freedom in the system upon decoding. 
Then corrections can be made accordingly. However, in 
contrast to the classical case, quantum errors occur in 
many different forms such as phase flips - and not just 
as bit flips. Furthermore, the quantum information must 
be preserved without ever measuring it, because measure- 
ment which obtains information about a quantum state 
inevitably disturbs it. There are excellent references on 
the theory of quantum error correction |^-|26|] . We limit 
the present discussion to detection codes only. 

For a code to detect errors, it suffices to choose the 
codeword space C such that all errors to be detected map 
C to its orthogonal complement. In this way, detection 
can be done unambiguously by a projection onto C with- 
out distinguishing individual codewords; hence without 
disturbing the encoded information. To make this con- 
crete, consider the code ]18[ 
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where the subscript L denotes logical states. An arbi- 
trary encoded qubit is given by 



\^)=a\{) L )+b\l L ) 
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After the four possible errors in Eq.()lO|), the possible 
outcomes are 
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with the /irsi order erroneous states \ipzi) and l^/z) 
orthogonal to the correct state \ipn). Therefore, it is 
possible to distinguish ( |l5| ) from ( |lrf ) and (0) by a pro- 
jective measurement during decoding, which is described 
next. 

The encoding and decoding can be performed as fol- 
lows. We start with an arbitrary input state and a ground 
state ancilla, represented as qubits A and B in the circuit 
in Fig. |. 
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FIG. 2. Circuit for encoding and decoding. Qubit A is the 
input qubit. H is the Hadamard transformation, and the sym- 
bol next to H is a controlled-not with the dot and circle being 
the control and target bits. I^i-s) are given by Eq s . ( j 2C| ) - (p2~|) . 
p4, ps are mixtures of the states in Eqs . (|lE|) - (|l 8|) and in 
Eqs.@-@. A phase error in either one of the qubits will 
be revealed by qubit B being in |1) after decoding, and in 
that case, qubit A will be rejected. 

To encode the input qubit, a Hadamard transforma- 
tion, H , is applied to the ancilla, followed by a controlled- 
not from the ancilla to qubit A (written as QVb a ). Let 
A and B be the first and second label. Then, the two 
operations have matrix representations 
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CN b , 



|^) = (a|0) + 6|l))|0) (20) 

|^) = ^(a|0)+6|l))(|0) + |l)) (21) 

1-03) = ^(a|0) + 6|1))|0) + (a|l) + 6|0))|1) (22) 

= i=(a(|00) + |ll)) + 6(|01) + |10)), (23) 



where Eq.(p3|) is the desired encoded state. 

The decoding operation is the inverse of the encoding 
operation (see Fig. ||) so as to recover the input (a|0) + 
6|1))|0) in the absence of errors. Phase errors lead to 
other decoded outputs. The possible decoded states are 
given by: 
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Note that the ancilla becomes |1) upon decoding if 
and only if a single phase error has occurred. Moreover, 
qubits A and B are in product states but they are clas- 
sically correlated. Therefore, syndrome can be read out 
by a projective measurement on B without measuring 
the encoded state. The decoding operation transforms 
the codeword space and its orthogonal complement to 
the subspaces spanned by |0) and |1) in qubit B, while 
all the encoded information, either with or without error, 
goes to qubit A. 

We illustrate the role of entanglement in the digitiza- 
tion and detection of errors as follows. Suppose the error 



is a phase shift on qubit A: |0) — > |0), |1) 
the encoded state becomes 
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The decoded state is now a superposition of the states 
given by Eqs.© and ©: 
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and the qubits transform as (Fig. 



Measurement of qubit B projects it to either |0) or |1). 
Because of entanglement, qubit A is projected to having 
no phase error, or a complete phase flip. 

We quantify the error correcting effect of coding using 
the discrete interpretation of the noise process, leaving a 
full discussion of the fidelity to Section III. Recall from 
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Eq.fllOD that the errors I ® J, I ®a z , a z ®I, and cr z ® a z 
occur with probabilities (1 —p a ) (1 — j>&) , (1 -p a , p a (1 - 
Pb) and PaPb respectively, and only in the first and the last 
cases will the output state be accepted. The probability 
of accepting the output state is (1 — Pa)0- ~ Pb) + PaPb 
whereas the probability of accepting the correct state is 
(1— p a )(l — pb). The conditional probability of a correct, 
accepted state is therefore 



(l-p„)(l- P6 ) 

(1 -p a )(l ~Pb) +PaPb 



1 - PaPb 



(31) 



for small p a , Pb- The code improves the conditional error 
probability to second order, as a result of screening out 
the first order erroneous states. 

We conclude this section with a discussion of some 
properties of the two bit code. First, the code also ap- 
plies to mixed input states since the code preserves all 
constituent pure states in the mixed input. Second, we 
show here that two qubits are the minimum required to 
encode one qubit and to detect any phase errors. Let C be 
the 2-dimensional codeword space and E be a non-trivial 
error to be detected. For phase damping, E is unitary 
and therefore EC is also 2-dimensional. Moreover, C and 
EC must be orthogonal if E is to be detected. Therefore 
the minimum dimension of the system is 4, which requires 
two qubits. However, using only two qubits implies other 
intrinsic limitations. First, the code can detect but can- 
not distinguish errors. Therefore, it cannot correct er- 
rors. Moreover, \ipiz) decodes to a correct state in spin 
A which is rejected. These affect the absolute fidelity 
(the overall probability of successful recovery) but not 
the conditional fidelity (the probability of successful re- 
covery if the state is accepted) . Second, the error a z cr z 
cannot be detected. This affects both fidelities but only 
in second order. To understand why these limitations are 
intrinsic, let {Ek} be the set of non-trivial errors to be 
detected. Since E^C has to be orthogonal to C for all 
k, and since C has a unique orthogonal complement of 
dimension 2 in a two bit code, it follows that all E^C are 
equal, and it is impossible to distinguish (and correct) 
the different errors. By the same token, for any distinct 
errors E^i and E}., Ek'EkC = C because they are both 
orthogonal to EyC, which has a unique 2-dimensional 
orthogonal complement. Therefore, a two bit code that 
detects single phase errors can never detect double errors. 
Finally, since a detection code cannot correct errors, it 
can only improve the conditional fidelity of the accepted 
states but not the absolute fidelity. Here, we only re- 
mark that the conditional fidelity is a better measure in 
our experiments due to the bulk system used to imple- 
ment the two bit code. A discussion of fidelity measures 
in our experiments and quantum error correction in bulk 
systems will be postponed until Sections |l| and VI . The 
system in which the two bit code is implemented will be 
described next. 



C. Bulk NMR Quantum Computation 

Nuclear spin systems are good candidates for quantum 
computers for many reasons. Nuclear spins can have long 
coherence times. Coupled operations involving multiple 
qubits are built in as coupling of spins within molecules. 
Complex sequences of operations can be programmed 
and carried out easily using modern spectrometers. How- 
ever, the signal from a single spin is so weak that de- 
tection is not feasible with current technology unless a 
bulk sample of identical spin systems is measured. These 
identical systems run the same quantum computation in 
classical parallelism. Computation can be performed at 
room temperature starting with mixed initial states by 
distilling the signal of the small excess population in the 
desired ground state. How NMR quantum computation 
can be done is described in detail in the following. 

The quantum system (hardware) In our two bit 

NMR system, |0) and |1) describe the ground and ex- 
cited states of the nuclear spin (the states aligned with 
and against an externally applied static magnetic field 
Bo in the +z direction). As in the previous section, we 
call the spins denoted by the first and second registers A 
and B. The reduced Hamiltonian for our system is well 
approximated by (% = 1) |p7|p8| (see also Fig. ||) 



UJ a UJb T 7T J 

U = — — <7 Z O I - —I ® a z + —a z 



(32) 



The first two terms on the right hand side of Eq. ( j32j ) are 
Zeeman splitting terms describing the free precession of 
spins A and B about the —z direction with frequencies 
w a /27r and uji,/2tt. The third term describes a spin-spin 
coupling of J Hz which is electron mediated. It is known 
as the J coupling. H. env represents coupling to the reser- 
voir, such as interactions with other nuclei, and higher 
order terms in the spin-spin coupling. 

Universal set of quantum logic gates A set of logic 
gates is universal if any operation can be approximated 
by some suitable sequence of gates chosen from the set. 
Depending on whether computation is fault-tolerant or 
not, the minimun requirements for universality are differ- 
ent. In the latter case, any coupled two-qubit operation 
together with the set of all single qubit transformations 
form a universal set of quantum gates p9j-|32|] . Both re- 
quirements are satisfied in NMR as follows. 

Single qubit rotations Spin flip transitions between 
the two energy eigenstates can be induced by pulsed ra- 
dio frequency (RF) magnetic fields. These fields, oriented 
in the iy-plane perpendicular to Bo, selectively address 
either A or B by oscillating at angular frequencies uj a 
or ujb- In the classical picture, an RF pulse along the 
axis r) rotates a spin about fj by an angle 6 proportional 
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to the product of the pulse duration and amplitude of 
the oscillating magnetic field. In the quantum picture, 
the rotation operator e~ % ^ C! ' ri rotates the Bloch vector 
(Eq.(BI)) likewise. Throughout the paper, we denote ro- 
tations of 7r/2 along the x and y axes for spins A and 
B by X a , Y a , Xb, and Yf, with respective matrix repre- 



sentations e" 



-i-J/iglcr, 



, and e 4 4 ' 



The rotations in the reverse directions are denoted by an 
additional "bar" above the symbols of the original rota- 
tions, such as X a . The angle of rotation is given explicitly 
when it differs from ir/2. This set of rotations generates 
the Lie group of all single qubit operations, SU(2). For 
example, the Hadamard transformation can be written 
as H = ie~ l ~ c,y e l ~5 am which can be implemented in two 
pulses. 

Coupled operations Quantum entanglement, essen- 
tial to quantum information processing, can be naturally 
created by the time evolution of the system. In the re- 
spective rotating frames of the spins (tracing the free 
precession of the uncoupled spins), only the J-coupling 
term, e^T- '*® ^ ^ is relevant in the time evolution. En- 
tanglement is created because the evolution depends on 
the state of both spins. A frequently used coupled "op- 
eration" is a time delay of jj, denoted by r, which cor- 
responds to the evolution e -' i -°"2® cr * . For instance, ap- 
pending t with the single qubit rotations e 1 ^"'® 1 and 
e if /igier z aDOU t _| axes f ^4 anc j ^ we implement 

the unitary operation 
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which is a cross-phase modulation between the two 
qubits. Together with the set of all single qubit trans- 
formations, x completes the requirement for universality. 
For instance, the controlled- not, CNba, mentioned in the 
previous section, can be written as (H®I) x (H<g)I), and 
can be implemented by concatenating the sequences for 
each constituent operation. It is also crucial that the free 
evolution which leads to creation of entanglement can be 
reversed by applying refocusing 7r-pulses, such that cre- 
ation of entanglement between qubits can be stopped. 
This technique will be described in detail later. 

Measurement The measured quantity in NMR exper- 
iments is the time varying voltage induced in a pick-up 
coil in the iy-plane: 

V(t) = -VoTr [ e- mt p(0)e mt 

X ( (ia x + o y )®I + I ® (ia x + a y ) ) ] . (34) 

The signal V(t), known as the free induction decay 
(FID), is recorded with a phase-sensitive detector. In 
Eq.(p3), the onset of acquisition of the FID is taken to 



be t = 0. If the density matrix p(0) has Pauli matrix 
decomposition p(0) = V, _ CijO~i <8> <Xj where 00,1,2,3 
are the identity matrix and o- x . VtZ respectively, then the 
spectrum of V(t) has four lines at frequencies + 4' 
!rr ~ f> Irr + 'f' !?r — 1' with corresponding integrated 
areas ("peak integrals") 
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Note that the expression cio — C13 (cio + C13) occuring in 
the high (low) frequency line of spin A is the coefficient 
of a x ® (|0)(0|) in p(0). Similarly, c 20 - c 23 (c 20 + 

C23) is the coefficient of a y ® |1)(1| (|0)(0|). They signify 
transitions between the states |0) <-> |1) for spin A when 
spin B is in |1) (|0)). Similar observations hold for the 
high and low lines of spin B. 
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FIG. 3. Energy diagram for the two-spin nuclear system. 
The transitions labelled ai ow , ahigh, bi ow and bmgh refer to 
transitions (|0) <-> |1))|0), (|0) <-> |1})|1), |0}(|0) <-> |1}) and 
|1}(|0) <-> |1}) respectively. 



Thermal States In bulk NMR quantum computation 
at room temperature, a pure initial state is not avail- 
able due to large thermal fluctuations (h~Lo a , hujb <SC kT). 
Instead, a convenient class of initial states arises from 
the thermal equilibrium states (thermal states). In the 
energy eigenbasis, the density matrix is diagonal with di- 
agonal entries proportional to the Boltzmann factors, in 
other words, pth = 2 e_M % wnere kT is the thermal en- 
ergy and Z is the partition function normalization factor. 
At room temperature, (-g^) ~ 10~ 6 , Z « Dim(p t /j) and 
Pth ~ (I — ■^)/Dim( / 9 t / l ) to first order. For most of the 
time, the identity term in the above expansion is omitted 
in the analysis for two reasons. First, it does not con- 
tribute to any signal in Eqs.(^5|)-(^8|). Physically, it rep- 
resents a completely random mixture which is isotropic 
and, by symmetry, has no net magnetization at any time. 
Second, the identity is invariant under a wide range of 
processes. Processes which satisfy £(I) — I are called 



G 



united. These include all unitary transformations and 
phase damping. Under unital processes, the evolution of 
the state is completely determined by the evolution of the 
deviation from identity (^J^), in which case the identity 
can be neglected. 

The unitality assumption is a good approximation in 
our system. The main cause of non-unitality is amplitude 
damping. In general, for a non-unital but trace preserv- 
ing process £, the observable evolution of the deviation 
can be understood as follows [B3[ . Rewriting 



p = vl + p s 



(39) 



where ps = p 
identity and v 



vl is the traceless deviation from the 
1/Dim(p). Then, 



£(p) = vl + p s 

p s = v(£(I) - I) + £( PS ) 



(40) 
(41) 



The observed evolution of the deviation is ps — > ps- The 
second term in Eq.(ftl|) comes from the evolution of ps 
when the identity is neglected, and the first term is the 
correction due to non-unitality. In our experiment, am- 
plitude damping is slow compared to all other time scales, 
therefore £(I) — I is small and can be treated as a small 
extra distortion of the state when making the unitality 
assumption. 

Temporal labelling One convenient method to cre- 
ate arbitrary initial deviations from the thermal mixture 
is temporal labelling The idea is to add up the 

results of a series of experiments that begin with differ- 
ent preparation pulses before the intended experiment, so 
as to cancel out the signals from the undesired compo- 
nents in the initial thermal mixture. Mathematically, let 
{Pfe} be the set of initial pulses and £(p) be the intended 
computation process. By linearity, ^2^£{PkPthPl) — 
£(Sfc PkPthPl)- Summing over the experimental results 
(on the left side) is equivalent to performing the exper- 
iment with initial state J2 k PkPthPl (on the right side). 
Temporal labelling assumes the repeatability of the ex- 
periments, which is true up to small fluctuations. 

Example As an example of the above theories, consider 



applying the pulse sequence in Fig. II C to the thermal 
state. 



CN 



Y 



X 



FIG. 4. Pulse sequence for CN. Time runs from left to 
right. 



The pulses are short compared to other relevant time 
scales. Therefore, other changes of the system during the 



pulses are ignored. The unitary operation implemented 
by the above sequence is given by 



CN = 







iff* 






" 1 - i 








1 



















1 + i 







1 - 


i 



(42) 



(43) 



similar to CN ba described in Section II C . 

The deviation density matrix of the thermal state is 
proportional to —7i. Neglecting the J-coupling term, 
which is much smaller than the Zeeman terms, 



Pth ~ -y o-* ® I + ~2 1 d 
= - Diag(w a + ui b , uj a - u bl 



(44) 

-uj a - uj b ) (45) 



where "Diag" indicates a diagonal matrix with the given 
elements. For simplicity, we omit the proportionality 
constant in Eq.(^), and rename the right hand side as 
pth- The sequence transforms pth to 



p cn = CNpthCN 



= - Diag(w a + Lu bl -u a 



2 



Y 



(46) 

LU b ,~UJ a +UJ b ,U! a - UJ b ) (47) 

(48) 



in which the populations of 1 01) and |11) are inter- 
changed. Therefore, CN and CN ba act in the same way 
on the thermal state. 

By inspection of Eqs.([44|) and d4S|), it can be seen that 
Pth and p cn have zero peak integrals given by Eqs.(|35|)- 
(|38|). To obtain more information about the states, a 
readout pulse X a can be applied to transform the two 
states to 



Pth = ~Y ct v <E)I+y Tt 



and 



Pen = rT ® 



u b 

— I < 
2 



(49) 



(50) 



Now, in p' th is a term a y ® I with coefficient C20 = — ^j- 
which contributes to two spectral lines at ^ ± ^ with 
equal and positive, real peak integrals. The readout pulse 
transforms the unobservable coefficient C30 in p t h to the 
observable — C20 in p' th , yielding information on the state 
before the readout pulse. Similarly, p' cn has a a y ®o~ z term 
with coefficient C23 — wn i cn gives rise to two spec- 
tral lines with real and opposite peak integrals (Fig. ||). 
All outputs in our experiments are peak integrals of this 
type carrying information on the decoded states. 
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FIG. 5. (a) Spectrum of A after a readout pulse on the 
thermal state, (b) Spectrum of A after QV and a readout 
pulse. 



III. TWO BIT CODE IN NMR 

We now describe how the two bit code experiment 
can be implemented in an ensemble of two-spin systems. 

are 



Modifications of the standard theories in Section [I B 
needed. These include methods for state preparation, de- 
signing encoding and decoding pulse sequences, methods 
to store the qubit with controllable phase damping, and 
finally methods to read out the decoded qubit. Fidelity 
measures for deviation density matrices are also defined. 

Spins A and B are designated to be the input and the 
ancilla qubits respectively. The output states of spin A 
are reconstructed from the peak integrals at frequencies 
w Q /27r± J/2. Fig. |^ schematically summarizes the major 
steps in the experiments, with details given in the text. 

Some notation is defined as follows. Initial states and 
input states refer to po and p\ in Fig. ||. The phrase 
"ideal case" refers to the scenario of having perfect logical 
operations throughout the experiments and pure phase 
damping during storage. 



control 
experiment 



coding i 
experiment 



Exptl 



Expt2 



Expt3 



Expt4 





Y a (8) 




PD 
















PD 










CN 


Y a (6) 


x a 










PD 










Y a (8) 




Ufa 




x a 










PD 












CN 


Y a (6) 


^enc 


u dec 




x a 




Pure state 
synthesis 


Input state 
preparation 


Encoding 


Phase 
damping 


Decoding 


Read out 



I'th 



FIG. 6. Schematic diagram for the two-bit code experi- 
ment. QV is used to prepare the initial state. Y a (6) is a 
variable angle rotation applied to prepare an arbitrary in- 
put state, which is then subject to phase damping (PD). In 
the coding experiment, encoding and decoding operations, 
Uenc and Udec, are performed before and after phase damp- 
ing, whereas in the control experiment, these operations are 
omitted. X a is used as a readout pulse on A to determine 
the output state ps in spin A. pi corresponds to \tpi) or pi in 
Fig. p. Details are described in the text. 



Initial state preparation It is necessary to initialize 
spin B to |0) before the experiment. This can be done 
with temporal labeling using two experiments: the first 
experiment starts with no additional pulses; the second 
experiment starts with QV (Fig. [II C| ). Therefore, the 
equivalent initial state is pth + Pen (all symbols are as 
defined in the example in Section ETO) : 







-L0 b 
-LU a + LO b 



+ (Jb I c3> a x . 







-UJ b 



(51) 
(52) 



The first term in Eq.(|52j) is the desired initial state. The 
second term cannot affect the observable of interest, the 
spectrum at w a /27r, because of the following. The iden- 
tity in A is invariant under the preparation pulse Y a (9). 
The input state is thus the identity, which has no co- 
herence to start with. Therefore, the output state after 
phase damping in both the control and the coding experi- 
ment is still the identity. This is non-trivial in the coding 
experiment. However, inspection of Eqs.(24)-(B7|) shows 
that spin A is changed at most by a phase in the coding 
experiment. While Eqs.(p4])-(|2"7|) applies only to the case 
when B starts in |0), the result can be generalized to any 
arbitrary diagonal density matrix in B (proof omitted) . 
It follows from Eqs.(|35|)-(|38|) that the second term is not 
observable in the output spectrum of A. 

In contrast, the input state in spin A can be a mixed 
state as given by the first term in Eq.(^2|), since the phase 
damping code is still applicable. Different input states 
can be prepared by rotations about the y-axis of different 
angles 6 € [0, tt] to span a semi-circle in the Bloch sphere 
in the xz-plane. Due to axisymmetry of phase damping 
(Eq.(|])), these states suffice to represent all states to test 
the code. 

We conclude with an alternative interpretation of the 
initial state preparation. Let the fractional populations 
of 1 00), |01), 1 10) and 1 11) be p 00 , p Q i, pio and pu in 
the thermal state. Then, the initial state after temporal 
labelling is 



(53) 



(poi +Pn) I ® + 2pio I® ]0> <0j 
+ 2(p 00 -pio) |00)(00|, 



where the identity term is not omitted, unlike Eq.(p2[). 
Temporal labelling serves to randomize spin A in the first 
term in Eq.(^3|) when B is |1). We have shown previously 
that the identity input state of A is preserved throughout 
both the coding and the control experiments in the ideal 
case. Consequently, only the last term in Eq. (p3|) con- 
tributes to any detectable signals in all the experiments, 
and we can consider the last term as the initial state. 
Having justified both pictures using the first term in 
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Eq.(p2[) and the last term in Eg. (|53|) as initial state, both 
will be used throughout rest of the paper. 

Encoding and decoding The original encoding and 
decoding operations are composed of the Hadamard 



transformation and CNba, as defined in Section II C. The 



actual sequences can be simplified and are shown in 
Fig. 0. 



Urn 



Uk = Y a x Y a X a Y 



FIG. 7. Pulse sequences to implement the encoder U enc 
and the decoder Udec- Time ru ns fro m left to right. X{Y) a n,\ 
and t are as defined in Section 



[IC 



The operator U enc can be found by multiplying the 
component operators in Fig. [j], giving 



u enc - v - 



1-10 

i i 

1-1 

1 i 



(54) 



The en code d states arc slightly different from those in 
section 



HE: 



|04 = i=(|00)+i|ll)) 
|U) = -L(i|01) + |10)) 



(55) 
(56) 



but the scheme is nonetheless equivalent to the original 
one. 

The decoding operation Udec is given by 



U de 



1 

V2 



1 -i 
-10 -i 
-i 1 
-i -1 



= uL 



(57) 



The poss ible decoded outputs are the same as in Sec- 
tion II B except for an overall sign in the single error 
cases. 

Storage The time delay between encoding and decod- 
ing corresponds to storage time of the quantum state. 
During this delay time, phase damping, amplitude damp- 
ing and J-coupling evolution occur simultaneously. How 
to make phase damping the dominant process during 
storage is explained as follows. 

First of all, the time constants of amplitude damping, 
Ti's, are much longer than those of phase damping, T2's. 
Storage times td are chosen to satisfy f<j < T% <C T\. This 
ensures that the effects of amplitude damping are small. 



The remaining two processes, phase damping and J- 
coupling, can be considered as independent and com- 
muting processes in between any two pulses since all the 
phase damping operators commute with the J-coupling 
evolution exp(— i a z ® a z nJtd/2). We choose Jt d to be 
even integers to approximate the identity evolution. As 
J is known with limited accuracy, we add refocusing 7r- 
pulses (3f| to spin B (about the y-axis) in the middle 
and at the end of the phase damping period to ensure 
trivial evolution under J-coupling. These pulses flip the 
z axis for B during the second half of the storage time so 
that evolution in the first half is always reversed by that 
in the second half. In this way, controllable amount of 
phase damping is achieved to good approximation. 

Control Experiment For each storage time td, input 
state and temporal labelling experiment, a control exper- 
iment is performed with the coding and decoding opera- 
tions omitted. Since phase damping and J-coupling can 
be considered as independent processes, and J-coupling 
is arranged to act trivially, the resulting states illustrate 
phase damping of spin A without coding. 

Output and readout For an input state prepared with 
Y a (0), the state after encoding, dephasing and decoding 
(/?5 in Fig. ^) is derived in Appendix ^ and is given by 
Eq.(A7) (from here onwards, uj a is omitted) 



p coded = [ CQS Q (1 _ ^ _ ^ + 2papb) ^ 

+ sin 9 (1 -p a ~p b ) a x ] <g) (I + a z )/2 
+ [ cos 6 (p a +p b - 2p a p b ) a z 

+ sin0 (-p a +Pb) a*] <8> (I-a g )/2. (58) 

In the control experiment, the corresponding output 
state is given by (Eq.jA^)) 

p control = [ CQS0 az + (1 _ 2p ^ sin0 Ux j % { j + az y 2 

(59) 

The initial state used in the derivation of Eq. (^9|) is the 
first term in Eq. (|52|) , and the encoding and decoding op- 
erations are as given by Eqs.(|54|) and (57). 

In the ideal case, the output state can be read out 
in a single spectrum. Recall that the coefficients of 
— <J y ® (J ± o z ) and — <j x ® (J ± a z ) are the real and 
imaginary parts of the low and the high frequency lines 
of A. Therefore, the coefficients of — a z ® (J ± a z ) and 
-a x ®(I±a z ) in pf ded and p f ntro1 can be read out 
as the real and imaginary parts of the low and the high 
frequency lines of A, if X a is applied before acquisition. 
This pulse transforms the ^-component of A to the y- 
component leaving the x-component unchanged, as de- 
scribed in Section II C. Note that only states with spin 



B being |0) (|1)) contribute to the low (high) frequency 
line. Therefore, in the coding experiments, the accepted 
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(rejected) states of A can be read out separately in the 
low (high) frequency line. There are no rejected states in 
the control experiments. 

The rest of the paper makes use of the following no- 
tation. "Output states" or "accepted states" refer to 
the reduced density matrices of A before the readout 
pulse, and are denoted by p c a oded = B (0\p c 5 oded \0) b and 
pcontroi = B (0|pg°"* ro '|0) B . Re j ecte H sta tes refer to 
s(l|/9§ oded |l)s from the coding experiments. 

The accepted and rejected states for a given input as 
calculated from Eq.(|58|) and Eq.([59|) are summarized in 
Table | 





z — component 


x — component 


input state 


cos 8 


sin 6 


coding expt.: 
accepted state 
rejected state 


(l-Pa-p6 + 2p a p(,) COS0 
(Pa+Pt-2p a p b ) COS0 


(1-Pa-Pb) sin6» 
(-Pa+Pb) sin6» 


control expt.: 
accepted state 
rejected state 


cos9 



(l-2p a ) sine 




1 

0.75 
0.5 
0.25 


-0.25 
-0.5 
-0.75 
-1 




0.25 0.5 0.75 
x-axis of Bloch sphere 




0.25 0.5 0.75 1 
x-axis of Bloch sphere 



FIG. 8. Predicted output states (a) with or (b) without 
coding. The arrows indicate the direction of time and the 
ellipses represent snapshots of the original surface of the Bloch 
sphere. 



TABLE I. Input and output states of spin A in the coding 
and the control experiments. 

The output states p c a oded and p™" irol , as predicted by 
Table |, are plotted in Fig. || in the xz-plane of the Bloch 
sphere of A. The north and south poles represent the 
Bloch vectors ±z (|0) and |1)). The time trajectories of 
various initial states are indicated by the arrows. The 
Bloch sphere is distorted to an ellipsoid after each stor- 
age time. We concentrate on the cross-section in one half 
of the iz-plane, and call the curve an "ellipse" for conve- 
nience. The storage times plotted have equal spacing and 
correspond to p a = 0,0.071,0.133,0.185,0.230,0.269. 
For each ellipse, pb is chosen to be the same as p a - The 
main experimental results will comprise of information of 
this type. 

Fidelity One can quantify how well the input states 
are preserved using various fidelity measures. In classical 
communication, the fidelity can be defined as the proba- 
bility of successful recovery of the input bit string in the 
worse case. In quantum information processing, when 
the input is a pure state, the above definition generalizes 
to the minimal overlap fidelity, 



T 



^(PoutPin) 



(60) 



We emphasize that Eq.(^0|) applies to pure input states 
only. For simplicity, fidelities for mixed input states will 



not be given here. The reason why Eq.(60) is sufficient 
for our purpose will become clear later. 

When pi n and p ou t are qubit states of unit trace with 
respective Bloch vectors fi n and r out , Eq. (pfj|) can be 
rewritten as 



T = xam? in -(l -I- r in ■ r out ) . 



(61) 



Recall from Eq.(|^) that, for phase damping, when fj„ = 

(^7^1/1^2)1 ^ out 



(e r X) e xt r y ,r z ). Therefore, 



Tin ' ^out 



—At/ 2 1 2\ . 2 

e (r x + rj + r z 



-Mrl 



(62) 
(63) 



where we have used the fact |fi„| 2 = 1 for pure states and 
p = (1 — e~ At )/2. The minimum in Eq.(|6l|) is attained 
for input states on the equatorial plane with r 2 x + r^ = 1 . 
Therefore 



,-At 



(64) 



With coding, the accepted state is (see Eqs.(24)-(E7J)) 



p coded = (X - p a )(\ - p b )p m + PaPbO- z p ln O- z . (65) 

If one considers the conditional fidelity in the accepted 
state, p ou t in Eq.(|60|) should be taken as the post mea- 
surement density matrix, 



Pout 



n coded 
ra 



n coded 
ra 



tT(p c a ° ded ) (1 - p a )(l - p b ) + p a p b 
(1 -PaPb)p in ~r PaPb&zPin&z ■ 



(66) 
(67) 



Note that the above expression is identical to the ex- 
pression for single qubit phase damping but with error 
probability p — p a pb- Therefore, coding changes the con- 
ditional error probability to second order, and the condi- 
tional fidelity is improved to Tc — 1 ~ PaPb- 

The amount of distortion can also be summarized by 
the ellipticities of the "ellipses" that result from phase 
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damping. The ellipticity e is defined to be the ratio of the 
major axis to the minor axis. Without coding, the major 
axis remains unchanged under phase damping, and the 



minor axis shrinks by a factor of e' 
Using Eq.pL 



therefore e 



T 



1 



(1 



(68) 



With coding, Tc is given by the same expression on the 
right hand side of Eq.(|68|). In the ideal case, the over- 
lap fidelity and the ellipticity have a one-to-one corre- 
spondence. In the presence of imperfections, the over- 
lap fidelity and the ellipticity, one being the minimum of 
the input-output overlap and the other being an average 
parameter of distortion, are more effective in reflecting 
different types of distortion. 

We now generalize to new definitions of fidelity for de- 
viation density matrices for the two bit code. In NMR, 
quantum information is encoded in the small deviation of 
the state from a completely random mixture. The prob- 
lem with the usual definitions of fidelity is that they do 
not change significantly even when the small deviation 
changes completely. This is true whether the fidelities 
are defined for pure or mixed input states. To overcome 
this problem, we introduce the strategy of identifying the 
initial excess population in 1 00} as the pure initial state 
so that usual definitions of fidelity for pure input states 
are applicable. This improves the sensitivity of the fi- 
delity measures and provides a closer connection to the 
pure state picture. 

The initial state in Eq. (|53|) can be rewritten as 



p = ap pure + (1 - a)p qmet , 
where a = 2(p o — Pio) = huja/^kT, and 



pP ure = I oo) (00 1 

quiet ^ ^ I" 



1 - a 



Ooi + Pu) I® |1)<1| 
+2 p w I®|0}<0|] 



(69) 

(70) 
(71) 



It has already been shown that pi ulet is irrelevant to the 
evolution and the measurement of pP ure when all pro- 
cesses are unital. Therefore pi met is neglected and the 
small signal resulting from the slow non-unital processes 
will be treated as extra distortion to the observable com- 
ponent. The input state prepared by Y a {9) can be written 
as 



Pin = apZ re + (l-a)p^. 



(72) 



For the state change pi„ — > £{pi n ), we consider the over- 
lap between p^ re and £(p^ re ) in place of the overlap be- 
tween pi n and £{pin)- This defines a new overlap fidelity 



X, 



, / pure c 1 / pu t'( \ \ 



similar to the pure state case. 



can be calculated from the experimental results in 
the following manner. The measured Bloch vector of A, 
r m , is proportional to that defined by s(0|£(/9^ re )|0) b- 
Due to limitations in the measurement process, this pro- 
portionality constant a is not known a priori. However, 
when 9 = in the control experiment, £(pf n ire ) = P^ re 
and f m = dfi n . Therefore, a = \r m \e=o can be deter- 
mined. In other words, |r m |6/=o is used to normalize all 
other measured output states before using the expression 
for Tg- 

The expression for T a can also be used for the con- 
ditional fidelity in the coding experiment if the post- 
measurement accepted output state is known. This re- 
quires tr(p™ ded ) = (1 — p a — p b + 2p a p b ) to be deter- 
mined for each storage time. The correct normalization 
is again given by the output at 9 = 0, which equals 
r m = tr(p c a ° ded )&r m . 

In summary, each ellipse obtained in the coding and 
the control experiment is normalized by the amplitude 
at 9 = 0: 



J-x = minf . — 



1 



|f m (0 = O)| 



(73) 



It is interesting to note that in contrast to the fidelity 
measure, the ellipticity measure naturally performs an 
equivalent normalization, and thus can be used for devi- 
ations without modifications. 

We now turn to the experimental results, beginning 
with a description of our apparatus. 



IV. APPARATUS AND EXPERIMENTAL 
PARAMETERS 

We performed our experiments on carbon-13 labeled 
sodium formate (CHOO-Na+) (Fig. HJ) at 15°C. The nu- 
clear spins of proton and carbon were used as input and 
ancilla respectively. Note that the system was heteronu- 
clear. The sodium formate sample was a 0.6 milliliter 
1.26 molar solution (8:1 molar ratio with anhydrous cal- 
cium chloride) in deuterated water |3(| . The sample was 
degassed and flame sealed in a thin walled, 5mm NMR 
sample tube. 



H— 13 C 







0"Na + 



FIG. 9. "C-labeled formate. The nuclear spins of the 
neighboring proton and carbon represent qubits A and B. 

The time constants of phase damping and amplitude 
damping are shown in Table [fi| The fact Ti < T 2 ensures 
that the effect of amplitude damping is small compared 
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to phase damping. The experimental conditions are cho- 
sen so that proton and carbon have almost equal T2's: 
the most theoretically interesting scenario in which iden- 
tical quantum systems are available for coding. Coding 
with qubits having very different T2's is described in Ap- 
pendix |b| 

The time constants of phase damping and amplitude 
damping are shown in Table |fi| The fact T x < T 2 ensures 
that the effect of amplitude damping is small compared 
to phase damping. The experimental conditions are cho- 
sen so that proton and carbon have almost equal T 2 's. 
This eliminates potential bias caused by having a long- 
lived ancilla when evaluating the effectiveness of coding. 
This also realizes a common assumption in coding theory 
that identical quantum systems are available for coding. 
Subsidiary experiments with qubits having very different 
T 2 's are described in Appendix 0. 





Ti 


T 2 


13 C 


9 s 
13.5 s 


0.65 s 
0.75 s 



TABLE II. Ti's and T 2 's for CaCl 2 -doped formate 
at 15° C, measured using standard inversion recovery and 
Carr-Purcell-Meiboom-Gill pulse sequences respectively. 

Phase damping arises from constant or low-frequency 
non-uniformities of the "static" magnetic field which ran- 
domize the phase evolution of the spins in the ensem- 
ble. Several processes contribute to this inhomogene- 
ity on microscopic or macroscopic scales. Which pro- 
cess dominates phase damping varies from system to sys- 
tem p^| . For instance, intermolecular magnetic dipole- 
dipole interaction dominates phase damping in a solu- 
tion of small molecules, whereas the modulation of direct 
electron-nuclear dipole-dipole interactions becomes more 
important if paramagnetic impurities are present in the 
solution. For molecules with quadrupolar nuclei (spin > 
1/2), modulation of the quadrupolar coupling dominates 
phase damping. Other mechanisms such as chemical shift 
anisotropy can also dominate phase damping in other 
circumstances. These microscopic field inhomogeneities 
have no net effects on the static field when averaged over 
time, but they result in irreversible phase randomization 
with parameters intrinsic to the sample. Another origin 
of inhomogeneity comes from the macroscopic applied 
static magnetic field. In contrast to the intrinsic pro- 
cesses, phase randomization due to this inhomogeneity 
can be reversed by applying refocusing pulses as long as 
diffusion of molecules is insignificant. 

Phase damping caused by the intrinsic irreversible pro- 
cesses alone has a time constant denoted by T 2 , while the 
combined process has a shorter time constant denoted 
by To. T 2 is measured by the Carr-Purcell-Meiboom- 
Gill ^5[ experiment using multiple refocusing pulses. T 2 * 
can be estimated from the line-width of the NMR spectral 



lines: during acquisition, the signal decays exponentially 
due to phase damping, resulting in Lorentzian spectral 
lines with line- width 1 /irT^ ■ 

In our experiment, X^'s for proton and carbon were 
estimated to within 0.05 s to be « 0.35 and 0.50 s. The 
storage times td were approximately 0, 62, 123, 185, 246, 
308 ms (n/J for n = 0,12,24,36,48,60). The max- 
imum storage time was 120 t, long compared to the 
clock cycle and was comparable to T 2 . The decay con- 
stant A, defined in Section [I A, was given by A = 1/T 2 . 
The error probabilities after a storage time of td were 
Pi = (1 - exp(— td/T£ i ))/2 for spins i = A,B. To recon- 
struct the ellipse for each storage time, 11 experiments 
were run with input states spanning a semi-circle in the 
iz-plane. Each input state was prepared by a Y a (9) pulse 
with 6 = mr/10 for n = 0, 1, • • • , 10. 

All experiments were performed on an Oxford Instru- 
ments superconducting magnet of 11.7 Tesla, giving pre- 
cession frequencies of wa/27t 500 MHz for proton and 
ws/27r w 125 MHz for carbon. A Varian UNITY Inova spec- 
trometer with a triple-resonance probe was used to send 
the pulsed RF fields to the sample and to measure the 
FID's. The RF pulses selectively rotated a particular 
spin by oscillating on resonance with it. The tt/2 pulse 
durations were calibrated, and they were typically 8 to 
14 /zs. To perform logical operations in the respective ro- 
tating frames of the spins, reference oscillators were used 
to keep track of the free precession of both spins, leaving 
only the J-coupling term of 195.0 Hz in the time evolu- 
tion. Each FID was recorded for 6.8 s (until the signal 
had faded completely). The thermal state was obtained 
after a relaxation time of 80 s (^> T-y 's) before each pulse 
sequence. 

Using the above apparatus and procedures, we per- 
formed the experiments as outlined in Section The 
experimental results are described in the next section. 
Tomography results following the evolution of the state 
are presented in Appendix O. 



V. RESULTS AND DISCUSSION 



A. Decoded Bloch spheres 



The output states, p c a oded and p% 



^control 



obtained as de- 
scribed in Sections III and IV , and the analysis that con- 
firms the correction effects of coding, are presented in 
this section. 

Fig. [l(] shows the accepted states in the xz-plane of 
the Bloch sphere of A. p% oded and p f ntro1 are plotted 
in Figs, [h] a and b. The ellipse for each storage time is 
obtained by a least-squares fit described later. 

The most important feature in Fig. [l^ is the reduction 
of the ellipticities of the ellipses due to coding, which rep- 
resents partial removal of the distortion caused by phase 
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FIG. 10. (a) Experimental data (circles) showing the out- 
put states from the coding experiment. Each ellipse (solid 
line) corresponds to one storage time and is obtained by a 
least-squares fit (Eq.([79|)) to the data. The storage times are 
n x 61.5ms for n = 0, 1, • • • , 5, and smaller ellipses correspond 
to increasing n. (b) Experimental data (crosses) and fitted 
ellipses (solid lines) for the control experiment. A replica of 
figure (a) is plotted in dotted lines for comparison. In both 
figures, uncertainties in the data are much smaller than the 
circles and crosses. 



damping - the signature of error correction. Coding is 
effective throughout the range of storage times tested. 

We quantify the correction effects due to coding us- 
ing the ellipticities. When deviations from the ideal case 
such as offsets of the angular positions of the points along 
the ellipses and attenuation of signal strength with in- 
creasing 9 exist, the minimum overlap fidelities and the 
ellipticities are no longer related by Eq. (|6"8|) . Since the 
ellipticity is an average measure of distortion which is 
less susceptible to scattering of individual data points, 
we first study the ellipticities. Moreover, since the devia- 
tions from the ideal case are small, we can still infer the 
fidelities from the ellipticities using Eq. ( |68| ) . A discussion 
of the discrepancies and the exact overlap fidelities will 
given later. 

Ellipticities In the ideal case, the ellipticity for each 
ellipse can be obtained experimentally as 



1(9 = 0) 



(74) 



where I denotes the intensity (amplitude square) of the 
peak integral. / is given by the x and ^-components of 
the output states as 



In the ideal case, 1(9) can be found from Table |: 
Icontroi(O) = 1 - 4(p a - pi) sin 2 9 (76) 

Icoded(9) = (1 - Pa - Pb + ZPaPb) 2 

-4p a p b (l -p a -p b + p a p b ) sin 2 9 , (77) 
and both are of the functional form 



hdeai(9) = A + B sin 2 9. 



(78) 



(75) 



Experimentally, the output Bloch vectors do not form 
perfect ellipses. We modify Eq.(|78|) to include signal 
strength attenuation with increasing 9 and constant off- 
sets in the angular positions: 

Iex P (0) = (A + Bsm 2 (9 + D))(l-C(9 + D)), (79) 

and perform non-linear least-squares fits of the experi- 
mental data to determine A,B,C and D. The fitted el- 
lipses plotted in Fig. [l(] follow from Eq. ( f79| ) and the fitted 
parameters. The ellipticities e are found using Eq.(]74|) by 
interpolating the intensities at 9 = and 9 — ^ . The el- 
lipticities are plotted in Fig. [n] a. The uncertainties of 
the fitted parameters originate from the uncertainties in 
the data, which are estimated to be « 1% for the am- 
plitude and 1.5 degrees for the phase in the measured 
peak integrals. These uncertainties are propagated nu- 
merically to the ellipticities as plotted in Fig. |ll] a. Ideal 
case predictions and simulation results are also plotted in 
Fig. [H] a. The simulation takes into account the major 
imperfection in the pulses and will be described later. 

Error correction The effectiveness of coding to cor- 
rect errors is evident when comparing the ellipticities 
from the coding and the control experiments (Fig. 11 a). 
Without coding, the ellipticity grows exponentially as 
e td/T£ a ( for fitted t0 be w o.4s). With coding, the 

growth is slowed down, with almost zero growth for small 
td- The suppression of linear growth of the elliptic- 
ity can be further quantified by weighted quadratic fits 
e = cq + citd + cit 2 d to the ellipticities. For the con- 
trol experiments, Co = 1.00 ± 0.01, c\ — 1.31 ± 0.21 and 
C3 = 8.8 ± 0.8 whereas for the the coding experiments, 
c = 1.10±0.02, ci = -0.24±0.29 and c 3 = 3.8±0.9 (t d 
in seconds). Therefore, the linear term "vanishes" due to 
coding. The small negative coefficient for the linear term 
originates from the scattering of the data point at zero 
storage time. 

To quantify the "cost of the noisy gates" caused by 
the imperfect pulses, we compare the ellipticities from the 
coding experiments and from the ideal case, the quadratic 
fits of which are respectively e c e ° d p f = 1.10-0.24^+3.80^ 
and e^taf = 1-00 + 0A5t d + 2.50^. The imperfections 
cause the ellipticity to increase by 0.1 at td = and this 
extra distortion decreases with td- We take advantage of 
the fact that the simulation results are close to the data 
points but are not as scattered to have a better estimate 
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FIG. 11. (a) Ellipticity and (b) inferred fidelity as a func- 
tion of the storage time in the coding and the control experi- 
ments. Error bars represent 95% confidence level. 



of this "cost of the noisy gates" . The simulation data can 
be fitted by e c °^ d = 1.06 + 0.32i d + 2.47^. Compared 
to the ideal case, the coding operations increase the el- 
lipticity by « 0.06 at td = 0, and this extra distortion 
remains almost constant for all td- 

The error probabilities as inferred from the ellipticitics 
p e = 1 — T t = |(1 — -) are plotted in Fig. pi] b as a 
function of storage time. 

Error correction is also manifested by expressing p e in 
the coding experiment as a function of the original p e 
in the control experiment, as plotted in Fig. [12[ The 
quadratic fit to the experimental results gives p c ° d ^ d = 



Co + cip + C2P where p stands for p e in the control ex- 
periment, c = 0.047 ± 0.008, ci = -0.05 ± 0.12 and 
C2 = 1.38 ± 0.40. Therefore, the expected improve- 
ment p — > PaPb is confirmed. Experimentally, the er- 
ror probabilities are larger than in the ideal case by 



at most 4.7% and these extra errors decrease with p. 
The quadratic fit to the simulation results (which is 
a good approximation of the experimental data) gives 
pf^ d = 0.032 - 0.032p + 1.783p 2 and differs from the 
ideal case by a constant amount of w 0.033 ± 0.003 for all 
p, which represents the cost of the noisy gates. In con- 
clusion, coding with noisy gates is still effective in our 
experiments, even though the noisy gates add a constant 
amount of distortion. 
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FIG. 12. Error probabilities in the coding experiments vs 
the corresponding values in the control experiments. Error 
bars represent 95% confidence level. The 45° line is plotted 
as a dotted line. 



B. Discrepancies 

While the data exhibit a clear correction effect, there 
are notable deviations from the ideal case. First, the 
ellipses with coding are smaller than their counterparts 
without coding. This is most obvious when the storage 
time is zero, in which case the coding and the control 
experiments should produce equal outputs. Second, the 
signal strength is attenuated with increasing 9 relative to 
ideal ellipses. Third, although the data points are well 
fitted by ellipses, their angular positions are not exactly 
as expected ( "^-offsets" ) . Finally, the spacings between 
the ellipses deviate from expectation. The causes of these 
discrepancies and their implications on error correction 
are discussed next. 

Gate imperfections: RF field inhomogeneity The 

major cause of experimental errors is RF field inhomo- 
geneity, which causes gate imperfections. This was de- 
termined by a series of experiments (details of which are 
not given here), and a thorough numerical simulation, 
as described below. The physical origin of the problem 
is as follows. The coil windings produce inhomogeneous 
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RF fields that randomize the angles of rotation among 
molecules. For a single rotation, the signal averaged over 
the ensemble decreases exponentially with the pulse du- 
ration to good approximation. A measure of the RF field 
inhomogeneity is the signal strength after a 7r/2 pulse. 
They are measured to be » 0.96 and 0.92 for proton and 
carbon respectively. In other words, a single 7r/2 pulse 
has an error of w 4-8%. 

RF field inhomogeneity affects our experiments in 
many ways. First, it attenuates the signal in both the 
coding and the control experiments, but the effects are 
much more severe in the coding experiments which have 
eight extra pulses. For instance, when the storage time 
is zero, the two experiments should have identical out- 
puts, but the ellipse in the coding experiment is actually 
5-15 % smaller. Second, for each ellipse, attenuation in- 
creases with 9 as the preparation pulse Y a (9) becomes 
longer. 

The effects of the RF field inhomogeneity are compli- 
cated, because the errors from different RF pulses are 
correlated, and the correlation depends on the tempo- 
ral separation between the pulses and the diffusion rate 
of the molecules. The correlation time of the RF field 
inhomogeneity is comparable to the experimental time 
scales. For this reason, predictions of the effects of RF 
field inhomogeneity are analytically intractable. 

Numerical simulations were performed to model the 
dominant effects of RF field inhomogeneity. We fol- 
lowed the evolution of the states assuming random RF 
field strengths drawn from Lorentzian distributions (also 
known as Cauchy distributions) with means and stan- 
dard deviations matching pulse calibration and attenua- 
tion for the 7r/2 pulses. All parameters in the simulation, 
including TJ's, were determined experimentally without 
introducing any free parameters. As the exact time cor- 
relation function for the errors was unknown, except for 
numerical evidence of a long correlation time, we assumed 
perfect correlation in the errors. The simulated ensemble 
output signal was obtained by Gaussian integration with 
numerical errors bounded to below 1.5%. The results 
were shown in Fig. 

Besides phase damping and error correction effects in 
the data, the simulations also reproduce extra signal at- 
tenuation in the coding experiments. The ellipticities 
obtained from simulations (see Fig. |j"T|) approximate the 
experimental values very well. Such agreement to exper- 
imental results is surprising in the absence of free pa- 
rameters in the model. Simulation results allow the dis- 
crepancy between the observed and the ideal ellipticities 
to be explained in terms of the RF field inhomogeneity 
and allow the "cost of coding" to be better estimated to 
be the constant 6 % increase in ellipticity or the w 3 % 
increase in error probabilities. 

The simulation results also predict increasing attenua- 
tion with 9. From the fitted parameter C (see Eq.([79[)), 
the amplitudes at 9 — n are w 4% weaker than the cor- 
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FIG. 13. Simulated output states, plotted similarly as in 
Fig. The simulation results are fitted by ellipses similar 
to those of the experimental data. 



responding values at 9 = in the simulations. Experi- 
mentally, this attenuation increases from « 8 to 15% (as 
storage time increases from to 308 ms) in the control 
experiments, and remains ss 8% in the coding experi- 
ments. Therefore, RF field inhomogeneity contributes to 
the attenuation but only partially. 

We conclude that RF field inhomogeneity as we have 
modeled explains the diminished signal strength in the 
coding experiments. The simulation quantifies the "cost 
of the noisy gates" . RF field inhomogeneity also explains 
part of the attenuation with increasing 9. We can also 
conclude that other discrepancies not predicted by the 
simulations are not caused by RF field inhomogeneity 
and these discrepancies are described next. 

Other discrepancies The simulation results show that 
RF field inhomogeneity does not explain why the attenu- 
ation at large 9 increases with storage time without cod- 
ing, and it does not explain the ^-offsets along the ellipses 
and the unexpected spacings between them. 

The increased attenuation with storage time at large 9 
can be caused by amplitude damping. A precise descrip- 
tion |37j of amplitude damping during storage is out of 
scope of this paper and we consider only the dominant 
effect predicted by a simple picture. Phenomenologically, 
the loss of energy to the lattice is described in the NMR 
literature by 



z(t) = z(oo) + (z(0) - z{oo)) e 



-t/T! 



(80) 



where z(oo) = 1 is the thermal equilibrium magnetiza- 
tion. As z(0) = ±1 at 9 = and n, we expect no changes 
at 9 = 0, but expect \z\ at 9 = ir to decrease by — 7 % 
for td ~ — 308 ms and T\ « 9 s for proton. Note that 
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refocusing does not affect spin A in the control exper- 
iments |Q but it swaps \0l) and halfway during 
storage, symmetrizing the amplitude damping effects in 
the coding experiments. Therefore, we expect increased 
attenuation with storage time in the control experiments 
only. This matches our observations that 1(9 = tt) de- 
creases from 8 to 16 % in the control experiments, and 
remains 8% in the coding experiments. Moreover, earlier 
data taken without refocusing (not presented) have the 
same trend of increased attenuation with storage time in 
both the coding and the control experiments. These are 
all in accord with the hypothesis that amplitude damping 
is causing the observed effect. 

The second unexplained discrepancy is that the out- 
put states span more than a semi-ellipse in the coding 
experiment but slightly less than a semi-ellipse in the 
control experiment. We are not aware of any quantum 
processes that can be a cause of it. It is notable that 
the output states and the fitted ellipticities can be used 
to infer the initial values of 9, and they are roughly pro- 
portional to the expected values for each ellipse. The 
proportionality constants are 5-8 % higher than unity in 
the coding experiment, and 0-1.6% lower in the control 
experiments. Moreover, similar effects are observed in 
many other experimental runs. Therefore, this is likely 
to be a systematic error. 

We have no convincing explanation for the anomalous 
spacings between the ellipses in the experiment. How- 
ever, from the fact that all the data points belonging to 
the same ellipse are well fitted by it, the anomalous spac- 
ings are unlikely to be caused by random fluctuations on 
the time scales of each ellipse-experiment. The effect of 
the anomalous spacings is reflected in the scattering of 
the ellipticities of the data, and the large uncertainties in 
the quadratic fits. 

While it is impossible to eliminate or to fully explain 
these imperfections, it is possible to show that the devi- 
ations cannot affect the conclusion that error correction 
is effective. 



as error correction. Therefore, the effects of error cor- 
rection can still be confirmed in the presence of all these 
small discrepancies. 



C. Overlap fidelity 

The two previous subsections dealing with the elliptic- 
ities provide an analysis of the global performance of the 
code. A stricter analysis is provided in this section using 
the overlap fidelities given by Eq. ((73) in Section [II . The 



minimal overlap reflects all defects and deviations hidden 
in the ellipticities as well as other distortions such as that 
caused by amplitude damping. 

All measurements are normalized using the amplitudes 
at 9 = as described in Section III . In the control exper- 
iments, the output states at 9 = are least affected by 
amplitude damping and RF field inhomogeneity. There- 
fore, the normalization can be done accurately. In the 
coding experiment, the signal attenuation at 9 = due 
to RF field inhomogeneity can lead to overestimated fi- 
delities. We determine the uncertainties due to RF field 
inhomogeneity by the following method. For each stor- 
age time, the amplitudes of the accepted and the rejected 
states at 9 = are summed. The sum is compared with 
the corresponding amplitude at 9 = in the control ex- 
periment to estimate the attenuation due to RF field in- 
homogeneity. The effects on the overlap fidelities are 
bounded to below 2%. The errors in the measured peak 
integrals are propagated to the fidelities which result in 
standard deviations no more than 0.7%. We apply sim- 
ilar procedures to the simulation results. The net error 
probabilities, given by 1 — T for the control and 1 — Tc 
for the coding experiments, are plotted in Fig. fL4l 
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Effects of the discrepancies We now consider the 
effects of the discrepancies on the ellipticities and the in- 
ferred fidelities in the experiments. First of all, radial 
attenuation of the signal due to RF field inhomogene- 
ity does not affect the ellipticities nor the inferred fideli- 
ties (taken as conditional fidelities). Second, different 
expressions for the "ellipticity" are not equivalent when 
the output states do not form perfect ellipses. However, 
they differ by no more than 7 and 3 % in the control 
and the coding experiments, ^-offsets along the ellipses 
are not reflected in the ellipticities, resulting in overesti- 
mated inferred fidelities. This is bounded by 3%. The 
scattering of the ellipticities due to anomalous spacings 
between the ellipses is averaged out with curve-fits to the 
data. The most crucial point is, none of these effects have 
a dependence on the storage time that can be mistaken 
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FIG. 14. Overlap fidelity as defined in Eq.(^). Points in- 
dicate experimental data and dashed lines indicate simulation 
results. Error bars represent 95% confidence level. 

The large difference in the rates of growth of error 
probabilities confirms the effectiveness of coding even 
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when a stricter measure of fidelity is used. 

VI. CONCLUSION 

We have demonstrated experimentally, in a bulk NMR 
system, that using a two bit phase damping detection 
code, the distortion of the accepted output states can 
be largely removed. These experimental results also pro- 
vided quantitative measures of the major imperfections 
in the system. The principle source of errors, RF field 
inhomogeneity, was studied and a numerical simulation 
was developed to model our data. Despite the imper- 
fections, a net amount of error correction was observed, 
when comparing cases with and without coding, and in- 
cluding gate errors in both cases. 

Our analysis also addresses several theoretical ques- 
tions in quantum error correction in bulk samples such as 
the fidelity measures of deviation density matrices. In the 
following, we conclude with some observations regarding 
quantum error correction in bulk systems, including syn- 
drome measurements, the equivalence between error cor- 
rection and logical labelling the applicability and 
advantages of detection codes and some issues in signal 
strength. 

Projective syndrome measurements traditionally em- 
ployed in the standard theory of quantum error correc- 
tion are impossible in ensemble quantum computation. 
Measurements via the acquisition of the FID do not re- 
duce individual quantum states and provide only "aver- 
age syndromes" . Moreover, the quantum states are de- 
stroyed after acquisition. However, the important point 
is that syndrome measurement is not necessary in error 
correction p9|Jl^ |. 

In each molecule, the syndrome bits carry the error 
syndrome for that particular molecule after decoding. 
These bits can either be used in a controlled-operation to 
correct the error |}9| , or in the case of a detection code, to 
"logically label" the correct and erroneous states. Con- 
versely, logical labelling to obtain effective pure states 
can be considered as error detection: unsuitable initial 
states are "detected" and are labelled as "bad" to start 
with. Both processes involve ejecting the entropy of the 
system to the ancilla bits. Error detection and logical 
labelling are therefore similar concepts. 

The distinction between correction and detection codes 
is blurred when using bulk systems. The objective of er- 
ror correction is to achieve reliable data transmission or 
information processing with high probability of success. 
When information is encoded in single systems, encoding 
with an error detection scheme will fail to provide reliable 
output for two different reasons: accepting an erroneous 
state or losing the state upon the detection of an error. 
Therefore, coding schemes capable of distinguishing and 
correcting errors are necessary to improve the probabil- 
ity of successful data processing. In contrast, in a bulk 



sample, a large initial redundancy exists upon prepara- 
tion, and the combined signal of all the accepted cases 
forms the output. Therefore, rejection of the erroneous 
cases results in a reduction of the signal strength in the 
improved accepted cases without necessarily causing a 
failure. Detection codes thus provide a tradeoff between 
probability of error-free computation and signal strength. 

As suggested by this analysis, and in concert with our 
experimental results, it thus makes sense to use detec- 
tion codes instead of correction codes in bulk quantum 
computing systems under certain circumstances. Fun- 
damentally, it is valuable to be able to interchange re- 
sources depending on their relative costs, This is illus- 
trated by the following simple example. Suppose a total 
pool of to qubits is available for transmission, and one 
just wants to correct for single phase flip errors of prob- 
ability p. Using a three bit code, one would obtain an 
aggregate signal strength of to/3, with fidelity 1 — 3p 2 , 
whereas with a two bit code, the accepted signal strength 
would be m(l — 2p)/2, with fidelity 1 — p 2 . Therefore, 
whenp < 1/6, the two-qubit code performs better in this 
model due to its higher rate. 

Another example relevant to bulk computation arises 
when the encoding and decoding circuits fail with prob- 
ability proportional to the number of elementary gates 
used. Although errors in consecutive gates can be made 
to cancel sometimes, this basic scenario is substantiated 
by our experiment, in which imperfect pulses contribute 
significantly to the net error. Assume now that we have 
n molecules, which are either two or three qubit systems. 
Let us compare the performance of the two and three bit 
codes, based on the strength of the correct output signal. 
Because the correction code requires at least three times 
as many operations as the detection code |i0| , the figures 
of merit obtained for the two schemes are n(l — 3p g ) and 
n(l — 2p)(l —p g ) respectively, where p g is the gate failure 
probability. In this model, the detection code performs 
better for p < p g /(l — p g ) due to the simplicity of the 
coding operations. 

A third example is the case of current state NMR quan- 
tum computation at room temperature, in which the in- 
trinsic signal strength decreases exponentially with the 
number of qubits [p|]s|]3"i|| . In this model, the initial sig- 
nal strength of an effective pure state of m qubits is ap- 
proximately of order 2 _m , and thus, for an ensemble of 
n molecules, the signal strengths of the outputs from 
the correction and detection codes are about n/8 and 
n(l — 2p)/4 respectively. According to this measure of 
performance, the detection code outperforms the correc- 
tion code for p < 1/4 (p < 0.27 in our experiments). 

If signal strength indeed decreases exponentially with 
to, then some interesting generalizations can be made. 
For arbitrary qubit errors, a i-error detection code has 
distance d > t + 1, while a terror correction code has 
distance d > 2t + 1 pa]. If one encodes k bits in the 
extra number of qubits used, I — k, satisfies the singleton 
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bound [g3|,|2|]4§), I - k > 2d - 2. Therefore, the out- 
put signal strengths for detection and correction codes 
would be approximately proportional to (1 — pf(p))/2 2t 
and l/2 4t , where f(p) is a polynomial. The detection 
code is thus always better asymptotically in this model 

This work illustrates how a careful study of dynamics 
in bulk quantum systems can provide a valuable oppor- 
tunity to demonstrate and test theories of quantum in- 
formation and computation. The development of tem- 
poral, spatial, logical, and related labeling techniques 
opens a window allowing information about the dynam- 
ics of single quantum systems to be extracted from bulk 
systems. Furthermore, by systematically developing an 
experimental toolbox of quantum circuits and quantum 
error correction and detection codes, experiments which 
test multiple particle quantum behavior become increas- 
ingly accessible. With improvements in the initial polar- 
ization in the systems and new labeling algorithms which 
do not incur exponential signal loss [Q, and with bet- 
ter methods to control the major source of error, the RF 
field inhomogeneity, we believe that further study of bulk 
quantum systems will complement the study of single 
quantum systems, provide new insights into the dynam- 
ical behavior of open quantum systems, and further the 
potential for quantum information processing. 
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APPENDIX A: MIXED STATE DESCRIPTION OF THE TWO BIT CODE 



Recall that the initial state after ancilla preparation is given by po = &z <& (7 + cr z ) /2 (see Eq. (|5^) , with to 
After Y a {9), the new density matrix is given by 

pi — (cos 9a z + sin 9a x ) ® (I + a z )/2 . 

Without coding, phase damping changes the density matrix to 

cos 9a z + (1 - 2p a ) sin 9a x <g> (I + a z )/2 . 

With coding, the encoding, phase damping, and decoding change the density matrix to P3, p4, and p§: 

P?, 



control 
P5 



„coded ^ CQS (g, az + ay (g, ax )/ 2 



+ shx6(a x ® I + 7® <Jy)/2 



p coded = CQs9 ( az g ^ + (1 _ 2p Q )(l _ 2p 6 )o-„ ® (T x )/2 

+ sin0((l - 2p a )a x ® 7 + (1 - 2p 6 )7 <g) CTj,)/2 
pf** = C os0a z (8 (7+ (1 - 2p«)(l - 2p 6 )a*)/2 
+ sm6>o- a ® ((1 - 2p a )7 + (1 - 2 Pb )a z )/2 

= cos 9a z ® (1 -p a -Pb + 2p a p b )(I + a z ) + (p a +p b - 2p a p b )(I - a z ) /2 
+ sin 6*^ (g) (1 -p a -p 6 )(7 + <r z ) + (-p a + p&)(7- cr z ) /2 

COS 9(1 -p a -Pb + 2p a p b )(J z + Sm0(l - p a -pb)(Jx ® (7 + £7a)/2 

+ cos (p a +p b - 2p a p b )a z + sin 0(-p a +Pb)&x ®{I-a z )/2. 
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Alternative understanding of the code 

The operation of the code can be further understood 
using the above mixed state description. Without coding, 
the input qubit is stored in the terms <j z ® I and a x ® 
/ which decay at different rates and distort the state 
asymmetrically. With coding, the input qubit is stored 
in the terms (a z <g> <r z + a y <g> a x ) and (a x ® I + I ® a y ) 
which decay at the same overall rate to first order. The 
resulting decoded state therefore shrinks radially up to 
first order. 

When subject to phase damping, the points on the 
Bloch sphere move inwards transversely without coding, 
distorting the shape axisymmetrically (Fig. |l5| ). In con- 
trast, the points move somewhat radially with coding, 
preserving the spherical shape better. 




0.5 

X axis Bloch sphere 

FIG. 15. Experimental data described in Section [y]. The 
solid lines, which are linear fits to the data points, indicate 
the flow directions of the points on the Bloch sphere subject 
to phase damping. In the no coding case, the flows are nearly 
planar, towards the z-axis, whereas in the encoded case, the 
flows are nearly radial, towards the center. 



APPENDIX B: THE CASE OF VERY 
DIFFERENT T 2 'S 

While the case of equal Ti 's is interesting from a theo- 
retical standpoint, different spins in a molecule typically 
have quite different Ta's. To study the two bit code in 
this regime, we performed experiments with carbon-13 
labelled chloroform dissolved in acetone ^5j,|). All pa- 
rameters were similar to the sodium formate sample, ex- 
cept for the relaxation time constants. 



In the chloroform experiment, Ti's were 16 s and 18.5 s 
and T2's were 7.5 s and 0.35 s for proton and carbon re- 
spectively. Separate experiments with the ancilla dephas- 
ing much slower or faster than the input were performed 
by interchanging the roles of proton and carbon. T 2 *'s 
and td's were as listed in The ellipticities are shown 
in Fig. 




J 1 



FIG. 16. Ellipticities obtained in the chloroform experi- 
ments, with (a) proton and (b) carbon as the ancilla. Carbon 
dephases much faster than proton. Error bars represent 95% 
confidence level. 

From Fig. [l6] a, it is apparent that coding almost re- 
moves the distortion entirely when a much better ancilla 
is available. The question is, is coding advantageous over 
storing in the good ancilla alone? Theoretically, coding 
is always advantageous because the error probability is 
always reduced from pi (i being the input spin) to p a Pb- 
Fig. [l6] b shows that experimentally, such improvement is 
marginal, because the advantage of coding is offset by the 
noise introduced. Therefore, when the T2's are very dif- 
ferent, the bottle neck is the dephasing of the bad qubit. 



APPENDIX C: TOMOGRAPHY RESULTS AT 
MAJOR STEPS 

Quantum state tomography H is a procedure to recon- 
struct the density matrix given a certain set of measure- 
ments. In 2-spin NMR systems, eight out of the fifteen 
coefficients, Cy, in the Pauli decomposition are obtain- 



able from the peak integrals, Eq.(35)-(p8|). The remain- 
ing 7 parameters can be obtained by repeating the mea- 
surement process with additional readout pulses before 
acquisition. These pulses permute the coefficients aj. A 
series of 9 experiments with different readout pulses is 
sufficient to reconstruct the complete deviation density 
matrix. 

We reconstructed the deviation density matrices in the 
coding experiments. The results for 9 = tt/2 and w 
123 ms are shown in Fig. The ideal matrices were 
calculated using equations derived in Appendix |A|. 
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FIG. 17. Experimental and ideal density matrices po, pi, pz, pi and p$. The basis is as indicated in the diagram. For 
each density matrix, the amplitudes (top) and the phases (bottom) for the corresponding entries are plotted. The amplitudes 
are shown in arbitrary units and the phases of entries with small amplitudes are omitted. The data are taken for a coding 
experiment with 9 — n/2 and td ~ 123 ms (24/ J). 
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